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1. Motivation 

The B- and D-meson leptonic decay constants fm, are important inputs into the comparison 
of theory with various experimental results in heavy-quark flavor physics. The decay rates for 
fully leptonic decays of B and D mesons arc proportional both to the decay constants fs,D and to 
the appropriate CKM matrix elements; extraction of the latter from experimental results therefore 
requires a precise determination of the decay constants. 

Further, the decay constants are also crucial for the prediction of rates for decay processes 
that are rare in the standard model, such as B s -+ /./ /j , now being measured to good precision 
by various LF1C experiments [1]. Rare processes such as this can be ideal testing grounds for the 
standard model, and are potentially quite sensitive to contributions from new physics [2]. Currently, 
uncertainty in the standard model prediction for this particular process is dominated by theoretical 
input, including fn ,, motivating a precision study of this quantity. 

Here we present an update of our ongoing analysis of the MILC asqtad ensembles [3, 4, 5] with 
the goal of extracting the decay constants fn.D- In addition to presenting the current state of the 
analysis, preliminary error estimates for the ratios /e//o and fitjfn, are presented. In combination 
with other precise numerical results for the D-meson decay constants (such as with HISQ charm 
quarks [8]), these may be used for improved determinations of the 5-meson decay constants. 

2. Simulation details 

We employ the MILC gauge-field ensembles with asqtad staggered fermions. The gauge en¬ 
sembles are generated with 2 + 1 flavors of dynamical sea quarks, and a tadpole-improved gauge 
action. Heavy quarks are included using the clover fermion action with the Fermilab heavy-quark 
interpretation. For further details on the lattice action and simulation parameters, see [6]. 


3. Analysis of correlation functions 


We compute heavy-light two-point correlation functions with both smeared and local sources, 
and both pseudoscalar and axial-vector interpolating operators. For a correlation function with 
source type i and sink type j, the expected functional form factorizes as 
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where N t is the temporal extent of the lattice and Nx denotes the number of excited states included. 
Making use of time-reversal symmetry, all correlators are considered only in the range 0 < t < 
N t / 2, with the results for t > N t /2 averaged in before the analysis is done. The decay constant is 
determined from the ground-state amplitude of the pseudoscalar-axial vector correlation function, 
obtained from a joint fit to four correlators: pseudoscalar-pseudoscalar with both sources local or 
smeared, and pseudoscalar-axial vector with both local and smeared pseudoscalar sources. 

For each two-point function, a corresponding fit range [finin+max] is chosen. We use t min /a = 4 
for all fits, and increase Nx until an acceptable correlated fit is obtained; this requires Nx — 3 on 
most ensembles, with Nx = 4 used for the ensembles with the largest N t . f max is held constant in 
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physical units, approximately 1.5 fm for fits to bottom-quark correlators, and 2.0 fm for charm. 
If this constraint gives t mSLX /a > 30, we set t max /a = 30 instead to avoid fit covariance matrices 
which are poorly determined by the available statistics; because the signal-to-noise ratio for the 
correlators drops quickly at large t, this is not expected to significantly affect the final fit results. 

With a large number of excited states included, several techniques are used to ensure the nu¬ 
merical stability of the fits. For all manifestly positive fit parameters, the logarithm of the parameter 
is used in the minimization in order to impose the positivity constraint on the fitter. In any multi¬ 
exponential fit, a potential ambiguity in the ordering of the energy states Eq,E\,E 2 ,... results in the 
existence of multiple equivalent minima in the space of all fit parameters. Rather than using the 
energies directly as fit parameters, we fit to the set 

Eq, logC Ei-Eq), log(E 2 -Ei),... (3.2) 

which imposes the explicit ordering Eq < E\ < £3 < ... on the fit function. 

Finally, we use empirical Bayesian constraints [7] for all correlator fits. Priors imposed on the 
excited states are relatively loose, and are used only for numerical stability. For the ground-state 
parameters, we employ a “two-stage" fitting procedure. In the first stage, a single-exponential fit 
is carried out restricted to data in a region corresponding roughly to an observed “plateau" in the 
effective mass m e ff(f); the plateau is taken from f max as given above, up to approximately 2.2 fm 
for B-meson fits, and 3.0 fm for D-meson fits. 

This first-stage fit is then used to set Bayesian prior constraints for the second-stage fit, which 
are carried out to the range [t m inTmax] as described above. The prior mean values for ground-state 
parameters are set to the first-stage best-fit values. Prior widths are set to three times the one-sigma 
classical errors on the parameters from the first stage. The results of this procedure are shown for 
a representative correlator are shown in Figure 1. 

4. Chiral-continuum extrapolation 

Using the results of the two-point fits and separate calculations of renormalization constants, 
we construct the renormalized decay amplitude 0# from the ground-state amplitude associated 
with the axial-vector sink: 


f 2 *«, = (4.1) 

where q and Q are the light and heavy quark making up heavy-light meson H q , respectively. This 
quantity is proportional to the decay constant times the square root of the meson mass, 0# = 
fn q ^M Hq . 

The chiral-continuum extrapolation for </>// is carried out separately for B and D mesons us¬ 
ing the framework of partially quenched heavy-light meson staggered chiral perturbation theory 
(PQHMs^PT) [3]. With the increased precision of this updated analysis, we find that to obtain 
acceptable chiral fits, we must include next-to-next-to-next-to-leading order (N 3 LO) analytic terms 
with cubic dependence on the light-quark masses. 

Our preliminary chiral-continuum fits are shown in Figure 2. We emphasize that these are 
preliminary, blinded results, and include an unknown multiplicative factor. The overall goodness 
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Figure 1: Results of the “two-stage" fit procedure described in the text, for correlators with a local pseu¬ 
doscalar source and either local pseudoscalar sink (red) or local axial-vector sink (green), shown as an 
effective mass. The band (black, dashed lines) shows the uncertainty used in the empirical Bayesian con¬ 
straint for the full fit, obtained from the “plateau" fit result for Eq as described in the text. The solid lines 
joining the effective mass points are reconstructed from the full fit, and show good agreement to low t. 


of fit is equal to ^ 2 /dof = 99.4/85 (p-value: 0.14) for the D-meson system, and ^ 2 /dof = 78.1/85 
(p-value: 0.69) for the D-meson system. 

5. Preliminary error budget for decay-constant ratios 

Of particular interest in this calculation are various ratios of heavy-light decay constants, since 
these quantities can be calculated with better precision than the individual decay constants due 
to the cancellation of various systematic effects. In our previous work [3], the ratios fnjfn and 
fsjI b were reported. In this updated analysis, we plan to also include results for the bottom-to- 
charm ratios /b s //d, and fn / /)> The combination of these ratios with more precise results for 
charm decay constants [8] will enable improved determinations of the bottom decay constants. An 
estimated preliminary error budget for each ratio is shown in Table 1; the procedure for estimating 
each uncertainty contribution is detailed below. 

We first consider systematic errors due to the renormalization. The lattice axial-vector current 
is renormalized with the combined factor 



In the heavy-flavor ratios fsj/ d s and fn/fn. the dependence on Z v 4 cancels, so only the ratios 

, /Z,/4 /Z,/4 and p, 4 /p» 4 contribute to the error. On the other hand, for the ratios fn / fn and 

y v bb cc n bq cq 

/bJ/b the factor Z V 4 > cancels; dependence on the ratio Z V 4 /Z v a remains, but because this quantity 
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Figure 2: Preliminary, blinded results for the chiral-continuum extrapolation of the decay ampitude r f 
for the D (upper panel) and B (lower panel) heavy-light mesons. Only full-QCD points (valence and sea 
light-quark masses equal) are shown in the plot, but the fit includes all available partially-quenched data. 


- which formally differs from 1 at order a s nr s a 2 - is found to have vary slowly with respect to the 
valence quark mass, so the overall contribution to the systematic error is expected to be negligible. 

For the ratios fn ,.//« and fnj fr>, effects of higher-order terms in the perturbative factor p A 4 

Qq 

enter only through the small variation of the perturbative renormalization factors with the valence 
light-quark mass; we retain the estimated error of 0.1% in the ratios as derived in [3], In the ratios 
fij/fi) and fn/fry - the light-quark valence mass dependence cancels, but dependence on the heavy- 
quark mass, which is expected to be significant, remains. We conservatively ignore correlations 
between the two p-factors appearing in these ratios, and estimate the projected systematic error 
by varying numerator and denominator by the value Pml x 0^> which p n |j lx = 0.1 and a s = a,y(2/a) 
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Table 1: Projected error budgets for the heavy-flavor decay-constant ratios. All errors are quoted as percent¬ 
ages. Error projections are preliminary, and may change when estimated from the full chiral fit in the final 
analysis. Entries marked with a dash show quantities that do not contribute due to cancellation in the given 
ratios. _ 


Source 

fD s /fD 

fsJfB 

fn/fo 

/b s //a 

Statistics 

0.2% 

0.4% 

0.6% 

0.5% 

Heavy-quark discretization 

0.7% 

0.6% 

0.8% 

0.8% 

Light-quark discretization 

0.2% 

0.1% 

0.3% 

0.3% 

Chiral extrapolation 

0.6% 

0.6% 

1.0% 

0.8% 

Heavy-quark tuning 

0.1% 

0.1% 

1.7% 

1.8% 

Z K 

0.0% 

0.0% 

— 

— 

Z,, 4 

V QQ 

— 

— 

0.2% 

0.2% 

Finite volume 

0.2% 

0.2% 

0.3% 

0.3% 

Higher-order p A 4 

0.1% 

0.1% 

4.1% 

4.1% 

Total projected error 

1.1% 

0.9% 

4.7% 

4.7% 


evaluated with a sa 0.045 fm (see [6] for notation.) 

Heavy-quark discretization errors are estimated via power-counting using the explicit formulas 
in Appendix A of [3]. Six such terms contribute at G[cr) and &(a s a)\ each term is multiplied by 
an unknown fit parameter, for which nominal values are taken from the preliminary joint chiral- 
continuum extrapolation shown above. 

Another source of uncertainty in the bottom-to-charm ratios is from the (mis)tuning of the 
heavy-quark masses; this error cancels almost entirely in the ratios of decay constants with the 
same heavy quark. To estimate the heavy-quark tuning error, we model the dependence of on 
the heavy-quark mass Mq by including the leading 1 /Mq dependence expected from heavy-quark 
effective theory, 

* Hq (M Q ) = *~-W (5 ‘ 2) 

Mq 

The parameters <f> „ and <j)f are estimated by comparing nominal values of </>// for bottom and charm 
from our previous complete analysis [3], Fractional uncertainties in the heavy-quark masses are 
estimated to be 1.8% for charm and 2.7% for bottom based on uncertainties in K reported in [6], 
and are then used with classical error-propagation formulas to estimate the uncertainty in the decay 
constants, assuming a correlation coefficient of 1 for the heavy-quark mass appearing in Mb and 
Mb s , and similarly for the charm mesons. 

Light-quark discretization errors are estimated to scale as fj{cx s tr), with improvement based 
on the reduction of the minimum available lattice spacing from a ~ 0.09 fm to a ~ 0.045 fm. Chiral 
extrapolation errors are estimated by taking the estimated errors from [3] and projecting them to 
reduce according to the lightest available light-quark mass in r\ units, r\m q . Statistical errors are 
projected to improve as y7V c f g (the number of sources used on all configurations except the coarsest 
is held constant at 8.) In terms of m n L, the largest volumes available are essentially the same as in 
the previous work [3], so the finite-volume error is estimated to be unchanged. 
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Adding the individual errors in quadrature, this leads to projected total errors on the percent 
level for the ratios fr>Jfo and fnfl/n- and total error of 4.7% for the bottom-to-charm ratios /g /fr> 
and fsj' /d s ■ The latter ratios, in combination with more precise determinations of the charm decay 
constants [8], could lead to values for /g and fly with precision comparable to other state-of-the-art 
calculations [9, 10, 11, 12] if the uncertainty due to higher-order perturbative renormalization can 
be reduced. 
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